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1. Introduction 

In the standard perturbative formulation of QCD the elementary excitations are 
coloured pointlike particles - gluons. At very small distances the gluons are asymptot- 
ically free while at distances comparable with the correlation length of the theory (the 
typical hadron radius), the effective interaction becomes so strong that the very concept 
of gluons looses its sense. The relevant degrees of freedom for the infrared regime of QCD 
are expected to be colourless extended objects (strings). The corresponding operators are 
the loop fields 

$(C) = ^ P exp[i j dx^(x)} (1.1) 

are defined in the space of all closed contours C. The letter P in front of the exponential 
stands for Dyson ordering along the contour. 

As it was pointed out by 't Hooft |J, the perturbative expansion in the iV-colour 
quantum chromodynamics has the same topological structure as the loop expansion in a 
string field theory with coupling constant l/iV0. One can imagine the world surface of the 
chromodynamical string as the result of condensation of dense planar diagrams. Then in 
the large iV limit the Wilson loop average W(C) = ($(C)) would be given by the sum over 
all planar surfaces S spanning the loop C, with some weight Q(S). Symbolically 

W(C) = ( L2 ) 
ds=e 

The perturbation theory does not tell much about the form of the Boltzmann factor 
Q(S) but it is known that the functional W{C) satisfies a loop equation involving a 
local deformation of the contour and having a local (vector) contact term. Therefore it is 
very plausible that there exists a string representation of the Wilson loop average of the 
type ( |1.2[ ) in which the weight of each surface depends on the geometry of the world sheet 
through a local action. 

The loop functional ( |1.2|) can be determined in principle from the following conditions. 
First, it does not change if a backtracking piece rr _1 is added to the contour C 

wicvv- 1 ) = w(C) (1.3) 

For such functionals ( called in Stokes type functionals) , one can define the so-called 
derivative S/Sa^ defined by 

6W(C) = {(dxpSxy - dxjxj !! W y. (1.4) 

1 Note that, unlike the usual string theories, the coupling constant is a dimensionless number 
and does not renormalize. 
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The equation of motion for the gauge field 



is transformed in the multicolour limit N — > oo to the following closed equation for the 
loop functional (|1.2| ) derived by Makeenko and Migdal @ 

It was shown in that the iteration of eq. (|1.6|) reproduces the series of planar diagrams, 
in its nonrenormalized form. Therefore this equation may help in principle to choose the 
right string theory for QCD. 

Unfortunately, eqs. ( |1.3| ) and (|1.6| ) are not of much practical use because they are 



written for regularized but not renormalized quantities. Therefore they imply a definite 
regularization of the corresponding string theory which may not be the simplest one @. 

The way out of this impass might be found by considering the lattice definition of 
the U(N) gauge theory The lattice formalism, in which the continuous Euclidean 
spacetime is replaced by a .D-dimensional simplicial complex, provides an unambiguous 
definition of both gauge theories and random surfaces. Once the lattice version of the 
string Ansatz (|1 . 2|) is knowm, one can try to formulate the world sheet dynamics of the 
QCD string. 

Some years ago, a mechanism of generation of strings within the strong coupling 
expansion in U( N) gauge theory was suggested by V. Kazakov [|| . The idea of |J , borrowed 
from the Stanley's solution || of the large- N vector model, was to reformulate the U(N) 
lattice gauge theory as a kind of Weingarten theory with additional interactions. The 
Weingarten model is a simple (but rather pathological) gauge theory which is equivalent 
to a model of noninteracting random surfaces in the large iV limit. In the U(N) lattice 
gauge theory one can impose the unitarity of the link variable by means of a Lagrange 
multiplier and replace the Haar measure on the U(N) group by a Gaussian measure over 
complex matrices. The construction suggested in || has been accomplished in refs. H, 
|J. In these papers it was shown that the strong coupling expansion of the Wilson U(oo) 
lattice gauge theory can be formulated in terms of a model of noninteracting surfaces on 
the lattice. 

Are these microscopic strings related to the strings expected to the strings expected to 
to arise in the continuum theory? Since the weak and the strong coupling phases in D = 4 



2 In the trivial case D = 2 there are no divergencies and the renormalized form of the loop 
equation ( |1.6| ) can be readily found and solved [H]. The solution was interpreted as a sum over 
minimal surfaces spanning the loop. 
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dimensions are separated by a first order phase transition, the mechanism of formation 
of strings might be very different in the two phases. We argued in ref. |T(J that the 
strong coupling lattice surfaces and the Feynman-'t Hooft planar diagrams correspond to 
two extreme cases of the same generalized model of random surfaces. The geometrical 
meaning of the phase transition is that in the weak coupling phase, windows (holes with 
free boundaries) may appear spontaneously on the world sheet. In the continuum limit 
the world sheet of the string is tore into a network of thin strips which can be interpreted 
as gluon propagators. In this way the weak coupling phase is achieved by relaxing the 
kinematical restrictions obeyed by the surfaces on the lattice. Therefore the strong coupling 
lattice string is perhaps not completely irrelevant to the continuum limit of the gauge 
theory. 



In this paper we develop further the approach outlined in |10[. We derive the string 



representation of the U(oo) gauge theory using a more traditional mathematical formalism 



known from the mean field analysis in gauge theories ||11|| . This formalism allows to 
construct the 1/N expansion in a systematic way. The leading term is given by the path 
integral of the free string, that is, by the sum over noninteracting random surfaces on the 
lattice. The intrinsic geometry of the world sheet of the string allows local singularities 
(punctures) with curvature 2nx integer. These singularities can occur at links or cells of the 
lattice. The Boltzmann weight of such surface is equal to a product of factors associated 
with the punctures and can have positive or negative sign i. These factors appear as the 
vacuum expectation values of special puncture operators located at the links and cells of 
the lattice. The world surface is also allowed to make folds but the total contribution of 
the surfaces with folds vanishes. 

The subleading terms in the 1/N expansion have double origin. First, there is the 
contribution of the surfaces with higher topology. Second, the fluctuations of the puncture 
operators, which are of order 1 /N 2 , lead to contact interactions of the surfaces. The vertices 
of these interactions can be imagined as microscopic tubes, trousers, ets. , connecting 
several punctures. 

We pay special attention to the case D = 2 where the equivalence between the gauge 
and string theories is established for a two-dimensional manifold M. with any topology 
discretized by a lattice C. It happens that in this simplest case this equivalence takes 
place up to the continuum limit which is the limit of an infinitely dense lattice L. 

Very recently, some interesting exact results have been obtained for the partition 
function of U(N) and SU(N) gauge theories defined on a compact two-dimensional surface 
of genus G, refs. [l^] [|14| [|15| [|16| ]17|]. An interpretation of these partition functions in 



terms of surfaces was suggested by Gross and Taylor [14|, [15]. The partition function 



3 This is the reason why the triviality theorem of Durhuus, Frohlich and Jonsson [12] is not 
applicable here. 
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can be written as a sum over oriented coverings of the space-time manifold of genus G 
and area A, where different sheets of a covering are glued together by means of branch 
points, microscopic handles and tubes, and what they called "omega points" . The "omega 
points" represent topologically branch points, handles and tubes but have no location on 
the surface. This representation looks particularly simple on the torus where there are no 
such singular points at all. On a surface with genus G there are |2 — 2G\ "omega points" 
and it is not clear how they can be incorporated in a local string Lagrangian. 

The comparison with the lattice string Ansatz allows to resolve the mystery about the 
"omega points" in the surface interpretation of the free energy. Indeed, in the continuum 
limit (dense covering lattice £), an "omega point" takes into account the contribution 
of a whole class of surfaces distinguished by the position of a localized singularity (e.g., 
branch point or microscopic tube). The sum of the Boltzmann weights of these surfaces is 
a topological invariant and is not sensitive to the choice of the covering lattice C. This is 
why the "omega point" has no fixed location in A4. 

The paper is organized as follows. 

In sect. 2 we give a review of the Stanley's solution || of the large- N vector model 
which may help the reader to understand the logic of our construction before dipping into 
the technical details. In sect. 3 we reformulate the functional integral for the gauge theory 
as an integral over complex matrix fields with flat measure and regular interaction. This 
can be achieved by introducing special puncture operators coupled to the moments of the 
matrix fields and playing the role of Lagrange multipliers. In sect. 4 we perform the 
integration over the matrix fields using the 't Hooft diagrammatic rules and interpret the 
result as the partition function of a gas of lattice surfaces interacting with the vacuum 
surfaces and among themselves through the puncture operators. In sect. 5 we study the 
case of D = 2 dimensions. We find explicitly the weights of the punctures and check the 
result on some simple examples of Wilson loops. Then we consider the two-loop correlator 
and the free energy on a compact lattice |17j where the contact interactions should be 
taken into account. In sect. 6 we discuss the form of the random surface Ansatz describing 
weak coupling phase and suggest a measure in the space of continuous surfaces for which 
the Wilson loop functional has well defined area derivative. In Appendix A we derive the 
microscopic loop equations for the puncture operators. 

2. The U(N) vector model 

In this section we remind how the vector field with global U(N) invariance can be 
represented in terms of statistical mechanics of random paths. We are using notations 
maximally close to those to be used later in the case of the U(N) gauge theory. 
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The fluctuating field in this model is an iV-component complex vector = 
(•u*, u^) defined at the points x of the .D-dimensional hypercubical lattice and hav- 
ing unit norm 

N 

<-u x = ^K| 2 = l (2.1) 

i=i 

We assume a nearest neighbour interaction /3u* ■ u y associated with the links i =< xy > 
of the lattice. Then the measure 

N N 

[du] =\\du i du i S(J2\<\ 2 - 1 ) =du*du5(u* -u-1) (2.2) 

i=l i=l 

and the interaction are invariant under global U(N) transformations 

u^Uu,; UeU(N) (2.3) 



The partition function 

Z = e NJ7 = I Y[[du x ] 11 e N P<- u y (2.4) 

x £=<xy> 

can be calculated order by order in the inverse coupling constant f3 by expanding the 
exponential as a series of monomials (strong coupling expansion). This leads to a diagram 
technique which in the limit iV — > oo involves tree-like clusters of loops and can be summed 
up explicitly in all orders. 

Alternatively one can transform, by introducing auxiliary fields, the integral over the 
original vector field into a duable Gaussian integral and find an effective field theory in 
which the number 1/N plays the role of a Planck constant. 

Introducing a new complex vector field h = (h % , i = 1,2, N) we represent the 
unimodular measure as a flat measure for the unconstrained complex vector field v = 
(v % , i = 1, N). We insert the exponential representation of the 5- function 

S(u j ,v j ) = (—) 2 dh l dh l ( h ' u,i (2 _ 5 ) 
2ixi J_ ioo 

in the definition ( |2.4j ) of the partition function for every component u J x of the vector field. 
Here hg and vg are unrestricted complex vectors on which the integration is performed 
with a flat measure 

N N 

dv = Y[ dv j dv j , dh = Y[ dh J dh J (2.6) 

3=1 3=1 
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Introfucing the effective potential F(h) by 

e NF W = J [du]e ^[«-.h+h-.u] (2 _ 7) 

we write the partition function ( |2.4[ ) in the form 

Z = <F T = III (dh x dv x e- N( - h *^+ h *<^ ]Je NF ^ J] e ^ v -< (2.8) 

The potential F(h) depends only on the radial degree of freedom 

a = h* ■ h (2.9) 

and can be expanded as 

oo 

a 



F ( a ) = £/»v (2 - 10) 

71=1 

The coefficients of the series are determined by the Ward identity 

(J2 N- 2 8/dh l d/dh* - l)e NF = (2.11) 

i 

or, in terms of the a-field, 

(—— + a——-l)e NF W-Q (2 12) 

In the limit N oo the solution is 

F(a) = VTT4^ - 1 - In 1 + v/1 + 4a (2.13) 

2 

Using the integral representation 

e NF(h*.h) = j daddte N[^. h - a)+ F( a )] (214) 

we finally write the partition function in a form allowing to perform exactly the integral 
over the vector fields 



Z = e 



NT 



' x (2.15) 

x e=<xx'> 
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Denoting by C the connectivity matrix of the D-dimensional hypercubical lattice 

D 

c = j2J2( ed/dx ' 1 +e ~ 9/9x ^ ( 2 - 16 ) 

we can write the result of the gaussian integration over the v, h fields as 

e NT = I l[da x da x e N ^ a ^- & ^ e N ^ = (e N ^) A (2.17) 

X 

where 

;F [6] = -^tr log(l -P&C) (2.18) 

X 

is the sum over vacuum loops of the vector fields. 

Formally the sum over vacuum loops can be obtained by expanding the exponential in 
the last line of (|2.15| ) as a sum of monomials and applying all possible Wick contractions 

(Kvyj) = 8 X J), (v*Jh yj ) = 5 x , y 5) (2.19) 

The vertex a x h* ■ h x is represented by a node at the point x and contributes a factor a x to 
the weight of the loop. The vertex /3v* v y is represented by a segment of the loop covering 
the link I =< xy > and contributes a factor (3. The weight of each loop is therefore a 
product of local factors (3 and a x associated with its nodes and links, correspondingly. 
Denoting by V any closed loop on the lattice, the r.h.s. of ( |2.18| ) can be written as 

^o[^]=^/3 #iinks n^ ( 2 - 2 °) 

r xer 

The action for the fields a, a is proportional to N and in the limit iV — » oo these fields 
freeze at their expectation values determined by the classical equations of motion 



or, equivalent ly, 



da 2a y ! 

a = — — (x\ -\x) 2.22 

da l l-aPC 



a + aa 2 = l, (x\ — -\x) = 1 (2.23) 

1-0&C 



The last equation has a simple geometrical interpretation in terms of a random walk on 
the lattice. It can be obtained also directly from the unimodularity condition u* • u x = 1. 
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Consider the two-point Greens function G xy = (u* • u y ). Repeating the same steps as for 
the partition function we find 



(v* • v„) = (xl-^ly) = J2&* nodes (3 #links (2.24) 
where the sum goes over all paths T xy between the points x and y. Therefore the second 



of the equations (|2.23|) is just the condition 

G xx = 1. (2.25) 

After diagonalysing the connectivity matrix we write it as an integral in the momentum 

space 

' (3 (2.26) 



^ 2 + E"=i[2sin(^/2)]2 

m 2 = — - 2D (2.27) 

a/3 

The continuum limit (infinite correlation length) is achieved when m 2 — > 0, i.e., 2D (3a = 1. 
When D < 2, the integral is infrared divergent when m — > and gives 

(3= f { d 9 % 2 l 2 = m D ~\ D < 2 (2.28) 
J (2-7r)^ p z + m 2 

In the marginal case D = 2 we find the asymptotic freedom law for the dependence 
of the renormalized coupling f3 = f3(m) as a function of the correlation length 1/m 

/?=-Lln^ (2.29) 

Above two dimensions the integral is convergent at m —>■ and the limit of zero mass 
is achieved at finite value of (3 = (3 C . In the weak coupling phase (3 > (3 C the equations of 
motion ( |2.26| ) have no solution which means that something in our construction is wrong. 
This contradiction is resolved if we admit that the vector fields acquire nonzero vacuum 
expectation values (v x ) = v c i, (h x ) = h c i determined by the saddle-point equations 

h c i=/3Cv c i = 2D/3v ch v cl = fih cl (2.30) 

which are solved by h c i = v c i = ( strong coupling £7(AQ-symmetric phase) or by m 2 = 
1 — 2D (3a = with h c i,v c i ^ ( weak coupling phase with spontaneously broken U(N) 
symmetry). Geometrically the weak coupling phase is characterized by possibility of the 
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the random walk to break into two pieces. Then the two-point Greens function is a sum 
of connected and disconnected parts 



G xy = v 2 cl + (x\ a \y) (2.31) 
1 — apC 



The unimodularity condition now reads 



1 f d D p 1 o P c 



Thus in the weak coupling phase there is a spontaneous magnetization v c i = y 1 — P c / P < 
1 and 2N — 1 Goldstone excitations with zero mass due to the spontaneous breaking of 
the U(N) symmetry. 

In D = 1 dimensions the vector model has no dynamical degrees of freedom and is 
analogous to the D = 2 gauge theory. In terms of the weak coupling parameter A related 
to by 2/P = sinh(A/2) we find 

a = tanh — , 8 = ; , (2.33) 

2' H 2sinhA/2 v ; 

and the two-point correlator is given by 

G xy = e -(V2)l*-i/l (2.34) 



This two-point correlator can be interpreted either as the sum over all random walks 
(with backtrackings) connecting the points x and y, or as the only minimal path (without 
backtrackings) connecting the two points, with a weight factor G\ = e~ x ^ 2 associated 
with each of its links. In the first interpretation the mean value of a is tuned so that it 
compensates completely the entropy due to backtrackings. In the case of the gauge theory 
the analog of the backtrackings are the folds of the world sheet of the string. We shall see 
that the D = 2 gauge theory can be formulated in terms of minimal nonfolding surfaces. 



3. The U(N) gauge theory 

In the Wilson-Polyakov formulation of the U(N) lattice theory [HJ the independent 
fluctuating variable is a unitary matrix = {(t^)], i, j = 1, iV} 

uju, = I (3.1) 
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associated with the oriented links £ =< xy > of the hypercubical latticed. The measure 
and the interaction are invariant with respect to local gauge transformations 



U< X y> ► U :E U <X y > U 2y (3-2) 

The amplitude U(C) associated with the parallel transport along the loop C = \l\t<2.---tn\ 
is given by the product of the link variables 

U(C) = U, 1 U, 2 ...U,„ (3.3) 

The action in the gauge theory is a functional on the loop fields (|3.3|) . The simplest 
nontrivial loops C are associated with the elementary squares of the hypercubical lattice 
which are also calles plaquettes. In what follows we will denote by U p the ordered product 
of the 4 link variables along the boundary of the plaquette p. Let S\ (U) be the one- 
plaquette action. It contains as a parameter the coupling constant A of the lattice gauge 
theory. The partition function is defined by the integral over all link variables 



/ [dU,]TTe s ^ (3-4) 

Je.es 



where 

JV 

[dV] = dUdutyutu - I) = n dUjdiufysfeiutyul - S>) (3.5) 

i,j=l k 

is the invariant measure on the group U(N). By T we denoted the free energy per degree 
of freedom. 

The one-plaquette action is subjected to the following two requirements. First, it 
should be a real function defined on the conjugacy classes of the group U(N). This means 
that it can be Fourier-expanded in the characters of the irreducible representations of the 
group U (N) 

e N2s ^= Xn(l)Xn(U)e N * S * (3.6) 

— oo<ni < . . . <njv <oo 

The character of the representation with signature n = {ni, un} depends on the eigen- 
values ui, Ujsf of the unitary matrix U as 

xa(u) = d w;_t;> . (3.7) 



4 We will use the notation l^ 1 =< yx > for the link obtained from £ =< xy > by reversing 
the orientation. More generally, if Ti and are two oriented lines, the product Til^ is defined, 
in case it exists, as the line obtained by identifying the end of I\ with the beginning of T2- 
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Second, we assume that the continuum limit is achieved when A — > and in this limit 

Sa(U)--^A 2 , U = e* A (3.8) 

The Fourier image of the action in this limit is proportional to eigenvalue of the second 
Cazimir operator C* 2 = — tr(U<9/<9U) 2 

\ A N 

St = -2M^(n) = -— (X)n? + $> - »i)) ( 3 - 9 ) 

i=l i<j 

We will assume that eq. ( |3.9| ) remains true also for finite values of A. The corresponding 
one-plaquette action 

&(U) = ^log J2 Xn(l)xn(U)e^/W = ± lo g (I\e^^ |U) 

— oo<ni < . . . <njv <oo 

(3-10) 

is known as /iea£ kernel action MM . The exponential of this action is a solution of the heat 



kernel equation on the group manifold 

[2^ + |(Ua/aU) 2 ]e^i U) =0 (3.11) 
The heat kernel action has the nice property of reproducing itself after group integration 



J [d\J}e N2S ^ (UiU)^ 2 ^ (UaU" 1 ) = e ^V 2 S Al+A2 (UxU a ). (-3^3) 

Therefore in two dimensions the heat kernel action has trivial scaling. 

Originally the lattice gauge theory was defined in terms of the action 

Sf(U) = ^!(U + ut) (3.13) 

known as the Wilson action. However this is not the best choice in our problem since its 
Fourier image is given by two different analytic expressions in the limit iV — ► oo, depending 
on the value of the coupling constant. For example, for the fundamental representation 
n / = [0,0,...,0,l] 

S "M^A, "1<1 < 3 ' 14 > 

The nonanalyticity at the point A = 1 known as Gross- Witten phase transition [^DJ is due 
to the fact that the integral for the inverse Fourier transform 



-i- / [cfU^W^u) (3.15) 



Xn(I) 

is saturated by the vicinity of a saddle point and the saddle-point solution can have two 
analytic branches. In the strong coupling domain A > 1 the eigenvalues of the unitary 
matrix U are distributed all along the unit circle, and in the weak coupling domain A < 1 
the density of the phases of the eigenvalues is supported by an interval [—(f), <fi] with <fi < tv. 
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3.1. The gauge theory in terms of nonrestricted complex matrices 



Now let us convert this model into a theory with flat measure and regular interaction 
using the Laplace transform of the measure in the fields as in the U(N) vector model. This 



has been actually done in the mean field analysis of lattice gauge theories [ fTTf . We insert 
the exponential representation of the 5-function 

6(U(,Vt) = (— ) 2 / dH)dH) e NBj(ui-vn+NHj(ui-vi) (316) 
2tvi J_ ioo 

in the definition of the partition function for every matrix element (Ui)l of the gauge fields. 
Here and are unrestricted complex matrices on which the integration is performed 
with a flat measure 

N N 

dv = n dv i ' dH = n dH j d Hj ( 3 - 17 ) 

For each oriented link I =< xy > only one pair of these fields is introduced and the 
convention 

V< yx > = V e -i = vj = V^ xy> , H <yx> = H e -i = hJ = H< xy> (3.18) 
is assumed. Denoting by 

e iV 2 F(H) = J [dU]e iVtr(Htu+Hut) ( g lg) 

the one-link mean-field integral we write the partition function ( |3.4|) in the form 
Z = e N ' T = [U (diI e dV £ e- Nt ^v e +ii e v]) 

r2 T?m.\ -i-r e N 2 S(V p ) 



(3.20) 



where V p denotes the ordered product of link variables along the boundary dp of the 
elementary square p 



ledp 



The measure is now flat and the integrand is regular. The constraint ( |3.1| ) is satisfied in 
the sense of mean values 

((Vjv, -I)...)=0 (3.22) 
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Assuming that the saddle point for the integral ( 3.20Q is the trivial field H = V = 



we can express the free energy as a sum of vacuum Feynman diagrams. The functional in- 
tegral can be calculated by expanding the exponentials and applying all Wick contractions 
between the H and V fields 

(muvht) = ((h})^)?) = ^ (3.23) 

3.2. The U(N) gauge theory a la Weingarten 

In order to give topological meaning of the 1/N expansion we need to expand the two 
potentials F and S\ in terms of the moments 

a n = £(HtH)», (3.24) 

Pn = |v*\ /3_ n = |vt- ; n = 1, 2, ... (3.25) 
of the matrix fields. The coefficients in the the expansions 



m = E £ (3-26) 

n=l ki,...,k n >l 



n=l k 1 ,...,k n ^O 

can be found from the corresponding Ward identities (Appendix A). 
The lowest-order coefficients are 



f[l] = 1 5 f[2] = "I /[1,1] = 1, /[3] = 2, /[4] = -5 

.— ) — A^sinh— 



S[i] = = e A/2 , s [2] = S[_ 2] = (cosh^) - iVsinh-^)e A , ( 3 . 2 8) 
— + 2(A^sinh— ' 



[i,i] = s [-i,-i] = [ - ^sinh— + 2(ATsinh— ) ]e , s^-i] 



The interpretation of the functional integral as a sum over surfaces is possible if 
the action is linear in the traces of the matrix fields. This is achieved by the integral 
transformation 

e A^(H) = I f[ {dand&n e N2 ^ * <hV-*£*]) c wM«] (3 2g) 

J n=l 



By means of another system of parameters we represent the exponential of the plaquette 
action in the form 

e N 2 5(v) = J f[ {d p nd p nd p_ nd p_ ne N 2 &(^ (3.30) 



n=l 
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We therefore introduce at each link I and at each elementary square p a set of auxiliary 
loop variables coupled to the moments of the matrix fields 

aM,&M; n = l,2,... (3.31) 

Pn{p)Jn(p); n = ±l,±2,... (3.32) 

We call them puncture operators by analogy with the matrix models for noncritical string 
theories. After that we are able to represent the integral ( |3.2U| ) as a theory of Weingarten 
type described by the partition function 

Z = e N ^ = ( / n^H^V,e- mr ( H J"^ +H ^ v J) 

1 . (3.33) 

f a n Wtr(H]H,)" Y[ e ^[/3 n (p)trV"+/3_„trvJ™]\ 

p;n>0 I a,(3 

where the Boltzmann weights of the surfaces are themselves quantum fields and the average 
with respect to them is defined by 

<•>**=/ II da4i)da n (£)e- N2§ -^^^ 

eeC -' n> ° (3.34) 
J! dUp)d^{p)e- N2 "^^ +S ^- 

The integration over the matrix fields in (J3.33[) will produce an effective action which 




will be interpreted in the next section as a sum over closed connected surfaces on the 
lattice. The leading term in this action is proportional to iV 2 . Therefore, in the limit 
N — > oo the integral over the scalar fields ( |3.31| ) and ( |3.32|) is saturated by a saddle point 
and these fields freeze at their vacuum expectation values 

(a n (£)) = &n, n (p)} = 0-n(p)} = Pn (3.35) 

The fluctuations of the scalar fields around their mean values are of order 1/N and their 
correlations will appear in the nonleading terms of the 1/N expansion. 

If we are interested only in the large iV limit, the fields ( |3.31| ) and ( |3.32| ) can be 
replaced with their vacuum expectation values even before the integration over the matrix 
fields. We therefore find the following expression for the free energy of the U(oo) gauge 
theory 



-TVtrCHjVf+Hfvj"; 



. x (3-36) 

"Q e f a„tr(Hj H ,r Y[ e V/3"(trV"+trvJ n ) 
?;n>0 p;n>0 

The parameters a n are determined by the condition of unitarity ( |3.22j ) and f3 n can be 
considered as coupling new constants defining the plaquette action. 
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4. Lattice strings (random surfaces) 

4-1. Feynman rules for the matrix fields 

By applying the rules of the Feynman- 't Hooft diagram expansion, the partition func- 
tion of the gauge model (|3.20|) can be reformulated as the partition function of a gas of 
lattice surfaces. An important feature of these surfaces is that they can have local singu- 
larities of the curvature not only at the sites, as it is the case in the original Weingarten 
model, but also at the 1-cells (links) and 2-cells (plaquettes) of the lattice. These surfaces 
are composed of elements with the topology of a disk, glued together by identifying pair- 
wise pieces of edges with opposite orientations along the links of the lattice (this is the 
meaning of the contractions (|3.23| )). There are two types of surface elements associated 
with the a and $ vertices. 

a) link-vertices 

The vertex Aa n (£)tr(HjH£) n will be represented graphically by a punctured disk 
with boundary going along the loop (££ _1 ) n . The H-matrices are associated with the 2n 
edges and the operator a n - with the puncture (fig. [j]). This surface element has zero area 
and Gaussian curvature 2n(n — 1) which we associate with the puncture. 

b) plaquette-vertices 

The vertex N /3 n (p)trVp corresponds to a punctured disk with area n covering n times 
the plaquette p and having as a boundary the loop (dp) n going n times around p (fig. ||). 
We will refer to this surface element as n-plaquette, or multiplaquette of order n. The 
V-matrices are associated with the edges of the boundary and the /3-operator - with the 
puncture. As before, the curvature 27r(?i — 1) will be associated with the puncture. 

The vertex A^/3_ n (p)trv| ?1 corresponds to the same surface but with the opposite 
orientation. Note that the /3 n -operator changes to /3_ n after changing the orientation. 
This means that the puncture is not just a point, but rather an oriented infinitesimal loop. 

Expanding the exponential in a sum of monomials, we can assign to each monomial 
a set of multiplaquettes of different orders covering the 2-cells (plaquettes) and a set of 
link elements associated with the 1-cells (links) of the lattice. Then we apply all possible 
contractions ( |3.23|) between V and H matrices. The geometrical meaning of these contrac- 
tions is identifying oppositely oriented edges of link and plaquette surface elements. The 
result is a collection of closed surfaces composed from multiplaquettes glued together by 
means of cyclic contractions (the link vertices). A piece of such surface is shown in fig. [| 
The weight of each surface is the product of the a and (3 operators associated with the 
punctures. 

Alternatively, one can imagine these surfaces as made only of n-plaquettes glued 
together along their half-edges. Such was the interpretation given in ref. ]1(J . 

The dependence on the number of colors N is, as in all N x N matrix models, through 
a factor N x where x is the total Euler characteristics. 
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The sum over all surface configurations is given by e F °\- ol ->p\ where ^[a, 0\ is the sum 
over all connected surfaces! 



Fo[a J] = NXiS) °M] ( 41 ) 



S:dS=0 

In the last formula x(<S) is the Euler characteristics of the surface and the Q-factor is the 
product of all a, /3-operators associated with the elementary disks from which the surface 
S is composed. The partition function of the gauge theory is obtained as the expectation 
value 

e T = ^f q [&M) (4.2) 

In this way the original gauge theory is formulated as a model of lattice random 
surfaces whose Boltzmann weights are products of local gauge-invariant puncture operators 
a n ,]3 n . It follows from ( |3.34[ ) that the interaction of these operators, mediated by the 
fields a n , fl n , is of order 1/N 2 . The vertices of this interaction are the coefficients in the 
expansions (|3.26| ) and ( |3.27| ) . They can be interpreted as microscopic surfaces with several 
punctures (fig. ^). 

The classical values of the puncture operators are determined by the saddle point 
equations 

I a „ = p = Vt H )«), -A. = p = i(V»); ia„ = j*, i/3„ = (4.3) 
n da n n n d(3 n n y n oa n n op n 

In the limit N —>■ oo the free energy is given by the sum of noninteracting random surfaces 
T = f Q (aJ) + J2F[a} + J2S x [(3]- ^ ±-a n a n - ^ ±p n /3 n (4.4) 

l p £,n>l P,n#0 

Applying the same arguments to the Wilson loop average W(C) = (|^V(C)) we find 
in the limiy iV — > oo a representation of W{C) as the sum over all planar surfaces bounded 
by the contour C. 

W(C)= Q W ( 4 - 5 ) 

S:dS=C 

The Boltzmann weight of a surface is the product of the mean values ( |3.35|) associated with 
the punctures. The O-factor of the surface S containing n-plaquettes, n = 1,2,..., 
and cyclic contractions of order n,n = l,2,..., is given by 



tyS] = Yl0n)^ P \*n)^ (4-6) 

n=l 



The trivial exponentiation is the main advantage of the present method compared with the 



character expansion used in [14], [15]. 
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4-2. Irreducible surfaces (a miracle) 

Unlike the random walk, the two-dimensional surface can exist in configurations with 
very uneven intrinsic geometry. The typical syngularity is a "neck" representing an inter- 
mediate closed string state of small length connecting a "baby universe" with the main 
body of the surface. The most singular configurations characterized by necks occupying 
a single link (fig. |5|) can be readily removed from the path integral by absorbing their 
contribution in the Boltzmann weights of the punctures. Namely, the sum over all re- 
ducible surfaces containing "baby universes" renormalizing the link vertices can be taken 
into account by adjusting the classical values a n the puncture operators located at links. 

In general, the classical fields a n , (3 n are complicated functions of the coupling constant 
A. Miraculously, the new Boltzmann weights become easily calculable and do not depend 
neither on the dimension D nor on the choice of the one plaquette action! This is why 
we will restrict the sum Tq over all connected surfaces to the sum Ti of irreducible ones 
which will render the random surface Ansatz much simpler than it its original version. 

Definition: 

A surface S immersed in a lattic C is reducible with respect to given 
link I G £ if it splits into two or more disconnected pieces after being 
cut along this link. A surface which is not reducible with respect to any 
I G C is is called irreducible. 
If we restrict the sum over surfaces to the irreducible ones, the free energy is obtained 
by the same formulae ( }4.3| ) and (O), with JF replaced by J-j. In what follows by sum 



over surfaces we will understand a sum over irreducible surfaces. This allows to simplify 
drastically the classical equations of motion ( fO| ) and solve them exactly. It turns out 
that, with this replacement, the expectation values of the puncture operators a n coincide 
with the coefficients /[ n ] in the expansion( |3.26| ). 

4-3. Evaluation of the link-vertices 



Instead of solving directly (|4.3| ) we would like to present a short-cut derivation using 
the hidden unitarity of the V-field. 

Consider the Wilson loop average for the closed contour C = ••) in the lattice C. 
It is given by eq. (O) in which the sum is restricted to irreducible surfaces. 



We will exploit the unitarity condition vvt = I, applied to the loop amplitude. It 
means that the Wilson average W(C) will not change if a backtracking piece ll~ x is added 
to the contour C 

W(Cir x ) = W(C) (4.7) 
The sum over surfaces spanning the loop Cll~ x can be divided into two pieces 

wicir 1 ) = w(c)w(ir l ) + w conn {cu- 1 ) (4.8) 
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The first term is the sum over all surfaces made of two disconnected parts spanning the 
loops C and il~ x . The second term contains the rest. The constraint ( f4.7| ) is satisfied 



for all loops if W(££~ x ) = 1 and W conn (C££~ x ) = 0. here is only one irreducible surface 
spanning the loop ££~ x ; it contains a single contraction ai(£) between £ and £~ x . Therefore 
&i = 1. 

Now consider a surface contributing to the second term Wi. The links £ and £~ x may 
be connected to the rest of the surface by means of the same link vertex or by two different 
link vertices. The condition that their total contribution is zero is 



n-l 
k=l 



+ ^a fc a n _ fc = 0, n = 2,3,... (4.9) 



Eq. ( fOl) is identical to the loop equation satisfied by the coefficients f n in the expansion 
( p.26| ), whose unique solution is given by the Catalan numbers (see Appendix A) 

*n = f H = (-) (n - 1) ^P^ (4-10) 

1 J n\[n — 1)1 

Since the evaluation of the link vertices is the most important point in the whole 
construction, we give a second derivation which is based only on the the sum over surfaces 
for the trivial Wilson loops 

wm-'r] = <^(v*vjr> = i, n=i,2,... (4.n) 

For each of these Wilson loops the sum over irreducible surfaces contains only finite number 
of terms, namely, the link- vertices contracting directly the edges of the loop ££~ x . For 
example, 

W(££- x ) = ai, W[(ir x ) 2 ] = 2a\ + a 2 , - (4.12) 
Introducing the generating functions 

oo 1 oo 

w{t) = Y,t n WW- x Y} = — v f(t) = l + J2t n (&n) n (4.13) 

n=0 n=l 

we easily find the relation |TtJ 



w(t) = f[tw 2 (t)} (4.14) 
which is solved by the generating function of the Catalan numbers 
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In this way, unlike the U(N) vector model, the classical values of the auxiliary fields 
do not depend on the dimension of the space-time. The expectation value a in the vector 
model is determined by the long wave excitations of the random walk. Here in the gauge 
theory, the expectation values of the puncture operators a n are determined in purely local 
way, as it is clear from their derivation, and coincide with the coefficients /[ n i in the 
expansion (|3.26|) of the potential F(a). One way to explain this difference is the local 



character of the U(N) invariance in the case of the gauge theory. 

On the contrary, the weights /3 n of the multiplaquettes will depend on the dimension 
as well as on the choice of the one-plaquette action. They coincide with the coefficients 
S[ n i in the expansion (|3.27| ) only in the trivial case D = 2. 



If we are interested only in the large iV limit, it is more convenient not to try to 
calculate them by solving the equations of motion but just to take them as independent 
coupling constants. 

We have found that after restricting the sum over surfaces to the irreducible ones, 
the expectation values of the puncture operators at the links are the same as the bare 
ones. That is, the renormalization of the bare link vertices /r n ] due to tadpole diagrams is 
neatly compensated by the contribution of the reducible surfaces. This is true also for the 
nonplanar vertices responsible for the contact interactions of surfaces. The equations for 
the renormalized link vertices obtained from the surface representation of the multiloop 
correlators (-^(H^hJ)^ 1 -^(H^hJ) 702 ...) = 1 are equivalent to the loop equations for the 
bare vertices considered in Appendix A. 

4-4- Loop equations 



In ref. [|10| we proved that the loop equations in the Wilson lattice gauge theory [^1 
are satisfied by the the sum over surfaces for the Wilson loop average. Let us only write 
here the general formula which is derived in the same fashion. For any link I £ T 

oo 

E E /3n[(W(C(^)-)-W(C(^)-^)] = E^(^)W(C^)[5(£,0-5(r\0] (4-16) 

n=lp:dp3t £'eC 

The sum on the l.h.s. goes over the 2(D — 1) plaquettes adjacent with the link I and the 
closed loops Civ , Cin in the contact term are obtained by cutting the links t and £' and 
reconnecting them in the other possible way. 

Eq. ( f!.16| ) is the lattice version of the loop equation ( |1.6| ) in the continuum theory. 
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5. The trivial D = 2 gauge theory as a nontrivial model of random surfaces 

The above construction can be easily checked in two dimensions where one can perform 
explicitly the sum over surfaces. The motion of the string is restricted by kinematics to 
the longitudinal vibrations and one can show that the contribution of the longitudinal 
modes of the string is zero, i.e. the surface maken no folds. The only singularities of 
a nonfolding surface are those associated with the punctures. A singularity containing 
curvature 2ii{n — 1) can be interpreted as branch points of order n. Therefore the Wilson 
loop average in two dimensions can be expressed as a sum over all surfaces with minimal 
area and eventually having branch points, bounded by the loop. This form of the random 
surface Ansatz in two dimensions was anticipated in and established in |TD| . 

5.1. The contribution of the surfaces with folds is zero 

The weights of the branch points, as they were determined by eq. (|4.9| ), provide a 
mechanism of suppressing the backtracking motion of the strings or, which is the same, 
world surfaces having folds. The surfaces with folds have both positive and negative weights 
and their total contribution to the string path integral is zero. 

We have checked that in many particular cases but the general proof is missing. An 
indirect proof might be constructed using the property ( |3.12| ) of the heat kernel action, 



which allows to delete parts of the lattice £ without changing the sum over surfaces. 

It is perhaps instructive to give one example. Consider the surface in fig. |6|. It covers 
three times the interior of the nonselfintersecting loop C and once the rest of the lattice. 
The change of the orientation of the world surface occurs along the loop. There are two 
specific points on the loop C at which the curvature has a conical singularity; they can be 
thought of as the points where the fold is created and annihilated. Each of these points 
can occur either at a site or at a link (in the last case it coincides with the branch point of 
order two of a link- vertex) . Let us evaluate the total contribution of all irreducible surfaces 
distinguished by the positions of the two singular points. Denoting by no and ni(= no) the 
number of sites and links along the loop C, and remembering that a singular point located 
at a link has to be taken with a weight 012 = /p] = — 1, we find that the contribution of 
these surfaces is proportional to 

n (n -l) 2 ni(ni-l) 

[ g + h\ n ^ + f [2 ] 2 1 ~ t ni + 2 -fe] n i] = ( 5 - 1 ) 

The second term on the left hand side contains the contribution of the reducible surfaces 
which must subtracted. A reducible surface arises when the two points are located at the 
extremities of the same link (m configurations) or when one of the points is a branch point 
and the other is located at one of the extremities of the same link (2ni configurations). 
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5.2. Evaluation of the plaquette-vertices 

Now we are ready to evaluate the expectation values of the fields $ n . For this pur- 
pose we consider the sum over surfaces for the one-plaquette Wilson loops W((dp) n ); n = 
±1, ±2, ... in the limit of infinite area A tot covered by the lattice C. In this limit the one-cell 
Wilson loops do not feel the rest of the lattice and coincide with these for the one-plaquette 
model. (This well known fact specific to two dimensions can be proved using the character 
expansion). These loops can be calculated from the loop equations (Appendix A); they 
coincide with the coefficients s n in the expansion ( |3.27| ). 

First let us observe that the only surfaces contributing to the sum are the surfaces 
with minimal area nX covering only the cell c. All other surfaces have folds (this is true 
only for an infinite lattice) and their contributions cancel. Therefore the sum over surfaces 
contains only finite number of terms and can be easily evaluated 

S[n] = W((dp) n ) = («i) 4n /3 n 

+ 4n(5i) nI - 2 (5 2 + («i) 2 ) faPn-k + - 

fc=i 

Applying the identities ( |4.9| ) one finds 

Pn = P-n = S{nY, U = 1,2,3, ... (5.3) 

A similar statement is true for the nonplanar plaquette vertices: amplitudes for the 
contact interactions of the puncture operators j3 n are given by the higher coefficients in 
the expansion ( [3.271) . 

5.3. String representation for the Wilson loops on the infinite plane. 

Now we are able to formulate the string Ansatz for the Wilson loop average as the 
sum over all irreducible surfaces without folds, spanning the loop C, and having branch 
points of all orders. It is convenient to write the O-factor in (|4.5|) as O = e _AArea O. Then 
the sum over surfaces (|4.5| ) takes the form 

W(C) = ^(S)e-5 AArea ( 5 ) (5.4) 
ds=e 

The factor fi(5) is a product of the weights of all branched points of the world sheet and 
the area of the surface is the sum of the areas of its elements (recall that the area of 
an n-plaquette is n). A branch point of order n (n = 1,2,3,...) associated with a /c-cell 
(k = 0, 1, 2 for sites, links, plaquettes correspondingly) is weighed by a factor where 

(o) 1 
w { = 1 



n 



(i) * < Un-i) (2n-2)! 



fn = (-) in - 1) 



n!(n-l)! (5.5) 



= s [n] e^ =y( \) ——r-{—^) m = (1 - ^— — ^-A + ...) 
H N \ m + l I m \ v ' v 2 

m=0 V 7 

(2) 

In the continuum limit A — > one can leave only the linear term in the weights u„ of the 
n-plaquettes. 
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5.4- Wilson loops on the infinite plane. Examples 

Let us illustrate how the string Ansatz ( |5.4| ) works by some simplest examples of 
planar loops (fig. [7]) which have been calculated previously using the Migdal-Makeenko 
loop equations 0. In order to simplify the notations, it is convenient to measure the 
area in units of 1/A. Then the dependence on A will completely desappear from the final 
expressions. 

a) Circle 

Assume that the loop encloses 2-cells and the total area is A = tt^A. The sum over 
surfaces contains a single nonfolding surface spanning the loop and 

W(C) =e~^ n2X = e~^ A (5.6) 

This is the famous area law for the Wilson loops in QCD 2 . 

b) Flower 

The nonfolding surfaces spanning the loop will cover the three petals of the flower 
(denoted by 1 in fig. |^) once, and the head (denoted by 2) - twice. Therefore the total area 
will be always A = A\ + 2A2, where Ai is the area enclosed by the domain i (i = 1, 2). 
Each of these surfaces will have a branch point located at some site, link or plaquette 
of the overlapping area 2 (fig. ||). Let us denote by no,ni,ri2 the numbers of sites, links, 
plaquettes belonging to the domain 2. Then the sum over the positions of the branch point 
reads 

^ Q = n k^[2] = ( n o - ni + (1 - X)n 2 ) = 1 - n 2 X (5.7) 

fc=0,l,2 

(Here we used the Euler relation no — n\ + n% = 1.) Therefore 

W(C) = (1 - A 2 )e-^ M+2M) (5.8) 

c) Pochhammer contour 

This time the surfaces spanning C may have different area. The three possible cases 
are shown in fig. [I[ There is one surface with area A\ + A^ 1 one surface with area A\ + A2 
and a whole class of surfaces with area A1+A2 + A2 which are distinguished by the position 
of a branch point of order 2. The two overlapping sheets are with opposite orientations 
and the branch point can occur only along the boundary of the overlapping region. Let 
no and n\ be the numbers of sites and links belonging to this boundary. By the Euler 
relation, no — n% = 1 — 2 = — 1 (the boundary has two cycles). Therefore the sum over the 
positions of the branch point contributes a factor J2k=o 1 UkU! [2] = n o — = —1 and the 
Wilson average reads 

W(C) = e~^ {Al+A2) + e ~^ Al+A ~ 2) - e" 5(^1+^2+^2) (5.9) 
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Finally, consider the sum over vacuum surfaces in the large iV limit. It follows from 
(fOp that the vacuum energy is equal to the sum over surfaces that, after being cut along 
all copies of a given link, split only into irreducible pieces. On the infinite plane, the 
only surfaces satisfying these conditions are made by identifying the boundaries of two 
multiplaquettes and the free energy is given by 

oo 2 

^=(#cells)^^ (5.10) 
^— ' n 

In the continuum limit A — ► 0, A tot = A^plaquettes = const, we find, up to an infinite 
constant, 

A 00 1 1 

^= j T i E^( 1 -- 2A + -) = ^tot (5.11) 

n=l 

5.5. Interactions 

The 1/N 2 corrections appeaar due to surfaces with higher topology as well as due to 
the contact interactions between the punctures. In two dimensions the interaction vertices 
are given by the coefficients in the expansions (|3.26|) and (|3.27| ). 



(AO 

[ni7 

n m located at a /c-cell (k = 1,2). Then, in two dimensions, 



Let u) , be the interaction amplitude for the puncture operators of orders 

[n 1 n 2 ...n m \ r f r 



^[1,1] = = s [i,i] ( 5 - 12 ) 

The simplest example where the interactions appear is the two-loop correlator 
W(Ci,C2)- Let us take the simplest case of non-selfintersecting loops (fig. [!(]). Then 
the answer depend on whether the two loops have equal or opposite orientations. 
/. Equally oriented contours 

(i) All cylindric surfaces bounded by the two contours have the same area, one cut 
and two branch points which can appear at the sites, links and cells of the overlapping 
domain 2. The sum over surfaces is a sum over the positions of these two branch points. 
The configurations leading to a reducible surface (i.e., when the two branch points can be 
connected by a link) are excluded from the sum. Let no, ni, ri2 be the number of sites, 
links and cells of this domain. Then the sum over the O-factors yields 

n (n - 1) -m (o) (o) , „ v (o) (i) 

2 U [2] U [2] + t n O n l - ^lj^[2] U {2] 

(0) (2) , 1 / -,\ (1) (1) 

+ noTtvU^u^ + -ni(m - l)ufa'w^ 

+ n!n 2 uj [2] u [2] H u [2] w [2] {o.ie>) 

= -x (no - ni + n 2 - 2An 2 )(n - ri\ + n 2 - 1) + ^A 2 n 2 (n 2 - 1) 
= ^A 2 n 2 (n 2 - 1) 
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(ii) The contact interactions due to the interaction of the puncture operators lead 
to a second term which is the sum over surfaces representing two discs connected with a 
microscopic tube. Their contribution is 

n 2U V> =n 2 (-X + ^X 2 ) (5.14) 

The sum over the two terms gives 

W(C X ,C 2 ) = e~^ +2A *\-A 2 + \a\) (5.15) 

II. Oppositely oriented contours 

(i) Besides the minimal cylindric surface with area A\ there is a whole class of cylindric 
surfaces with area A\ + A 2 (see ref. InJ). The contribution of these surfaces is ( [flOf , 
Appendix B) 

e~* Al + (n 2 - l)e~^ Al+A2 (5.16) 

(ii) The contact interactions produce a term 

-n 2 e~^ M+M (5.17) 
W(d,C 2 ) = e -3(^i-^) - e -h^+A 2 ) (51g) 
so that the total sum depends only on the two areas 

W(C) = - e -? Al+A2 + e~? Al (5.19) 

5.6. Compact lattices 

For a compact lattice £ with the topology of a sphere and area A tot the sum over 
surfaces becomes much more involved. It will contain a sum over spherical surfaces repre- 
senting multiple branched coverings of C. There are infinitely many surfaces contributing 
to the elementary Wilson loops. For example, 

W 1 = e~^ x + e~^ (Atot) + ... (5.20) 

Therefore, we can either modify the weights of the irreducible surfaces, or use the same 
weights as for the infinite lattice but allow contact interactions due to microscopic tubes, 
trousers, etc., connecting two, three, etc. multiplaquettes. The amplitudes of these contact 
interactions are given by the coefficients in the expansion ( p.27|) . 

We are not going to discuss this case in details since there is nothing conceptually 
new. As a single example we will calculate the two leading orders of the free energy due 
to the surfaces defining simple and double coverings of C 
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Let rii be the number of i-cells of the lattice C (i = 0, 1, 2). The free energy should 
depend only on the total area A tot = r^A. The two simple coverings (with the two pos- 
sible orientations) contribute the term 2e~^ Atot . There are two kinds of double coverings 
depending on the orientations of the two layers. 
Equal orientations 

The sum over cylindric surfaces with a cut is given by the same formula as for the two-loop 
correlator, with the only difference that here the Euler relation gives no — n\ + n 2 = 2. 
The result is (1 — Xn 2 + ^X 2 )e~ XAtot . We have to add the term with contact interaction 
due to the vertex wl 2 ^ = s^i] = (—A + ^A 2 ) connecting two equally oriented plaquettes: 
77,2(1 - X)e~ Atot . 
Opposite orientations 

The sum over irreducible surfaces consisting of two sheets with opposite orientations and 



connected along a tree-like configuration of links, was calculated in Appendix B of [10| 



and equals (p — 2)e~ Atot . We have to add to it the term with contact interaction due 
to the vertex u>^ ^ = /[i,-i] = —1 connecting two plaquettes with opposite orientations: 
— n 2 e~ Atot . 

Collecting all terms, we find 

T = 2e-^ + [-1 - 2A to t + ^ 2 ot]e" Atot + - (5.21) 

An interesting phenomenon on compact lattices was observed recently by Douglas and 
Kazakov ||17|| . The sum over surfaces for the free energy on the sphere is convergent only 



for sufficiently large total area 

A tot > 7T 2 (5.22) 

At the critical point A to t = tt 2 the entropy of the branched points overwhelms the energy 
and the effective string tension vanishes. The phenomenon resembles to what happens 
with the compactified bosonic string at the Hagedorn radius. 

The singularity at A = tt 2 can be observed already in the behaviour of the multiple 
Wilson loops surrounding a contour of area A. From the explicit expression given in the 
Appendix (with A replaced by A) one finds the asymptotics at n — > 00 

W n (A) « -^=(nVA)- 3 / 2 cos(2nv / I - 3tt/4) (5.23) 
V 71 " 

The value of A for which the oscillations disappear is exactly the critical area A c = n 2 . 
This is not a surprise since on a compact lattice with area A, the free energy contains a 
term J2 n Wn(A) (compare with eq. ( |5.1Q )). 

Finally, let make a remark concerning the case of a toroidal manifold M. . It was nicely 
demonstrated by Gross and Taylor that the free energy of the gauge theory on the torus 
can be interpreted as the sum over all connected nonfolding oriented surfaces wrapping the 



25 



torus, with right combinatorial factors. If we add to this sum the sum over microscopic 
surfaces ( |5.10|) , the result will be exactly the logarithm of the ^-function 



torus = — Atot - V log(l - e~ nA ™) (5.24) 



12 Wl ^ 

n=l 



6. Strong versus weak coupling 

6.1. The lattice string in the weak coupling phase 

The assumption that H c i = V c i = is a local minimum of the free energy is justified 
in the strong coupling phase of the gauge theory. In the case of Wilson action ( |3.13|) it 
is known that the weak coupling phase the matrix fields will develop vacuum expectation 



values |l 1J . The classical fields form an orbit of the gauge group 

V <xy> = U^VdU- 1 , H <xy> = U^HeiU- 1 (6.1) 
and are solutions of the saddle point equations 

a^(Vd) _ ^(Hci) (62) 



In the case of the Wilson action ( 3.13 ) these equations have only trivial solution H c i = 
V c i = for A > Xb and a nontrivial solution at A < Xb- It is quite possible the above 
is true for the general one-plaquette action. In the mean field approximation the critical 
coupling A c > As is determined by the assumption that the system chooses the absolute 
maximum of the classical free energy. 

Note that the saddle point equations does not fix the local gauge transformation 
{U x ;x G £}. The entropy of the gauge degenerate mean field ( |6.1| ) must be taken into 
account. In the continuum limit it produces zero modes among the oscillations around 
( |5~T| ) which has to be elliminated by a gauge fixing procedure. 

The diagram technique for the weak coupling phase is obtained by expanding the 
effective action S\(V) + F(H) — |VH around the classical solution ( |6.1|) 

H, = H cl + 5H,; V, = V cl + 5V, (6.3) 

The propagator will remain the same while new vertices will be generated. They will 
correspond to new surface elements having one or more free edges associated with the 
classical field H c i or V c i while the rest of the edges are associated with the quantum 
fluctuations 5V and 5H. After gluing pairwise all edges labeled by the fluctuations, the 
free edges will form closed loops. Geometrically a planar diagram composed by these 
elements represents a surface with free boundaries (windows) corresponding to these loops 
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(fig. |TT|). There are two kinds of windows: those made out of H c i edges and those made 
of V c i edges along their boundary. 

The gauge transformation in the classical solution ( |6.1|) can be absorbed in the inte- 
gration measure of the fluctuations. Therefore the Boltzmann weights of these generalized 
surfaces are gauge invariant. 

Let us also mention that in the weak coupling phase the potential for the H field is 
given by a different expression which is not analytically connected with the strong coupling 
solution ( |3.26| ) (see Appendix A). 

If the argument is proportional to the unit matrix, H = HI, then the weak coupling 
solution gives 

F(HI) =2H + const - I lnff + + -Lj + ... (6.4) 

We have seen that in the weak coupling phase characterized by a nontrivial classical 
solution, windows will appear spontaneously on the world sheet of the string. This situation 
resembles what happens in the weak coupling phase of the vector model in which the U(N) 
symmetry is spontaneously broken and which is described by random walks allowed to have 
free ends. 

Whether the windows destroy or not the world sheet in the continuum limit is a 
dynamical question. Believing in the confinement, we expect that in D < 4 dimensions 
the windows are still not sufficiently large for that, and the world sheet will have in the 
continuum limit the structure of a dense network of thin strips separating the windows. 
These strips will correspond to the gluon propagators in the standard Feynman rules. At 
the critical dimension D = 4 the effective string tension of the string (with windows on 
the world sheet) should scale with the coupling A according to the asymptotic freedom 
law, just as the mass in the vector model does in D = 2 dimensions. In D > 4 dimensions 
the world sheet of the string is eaten by one or several large windows and the Wilson loop 
behaves as exponential of the length of the contour. 

The above qualitative picture was based on the known mean-field solution for the 
Wilson action. Whether it is true in geleral, can be decided by studying the mean field 
problem for a general one-plaquette action. 

Such kind of world surface with windows is too complicated as an object on the lattice. 
However, if the string hypothesis is true, its critical fluctuations can be described in terms 
of an embedded surface with continuous world sheet and measure defined by a local action. 
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6.2. About the path integral for the continuum string 

Our conjecture is that the embedded surface is defined by two fields, the coordinate 
x n(£)i A* = 1 5 •••) D, and a field taking its values in the Grassmann manifold (jd,2 of 
oriented two-dimensional planes in Hp. The field counts for the additional degrees 
of freedom related to the tangent planes along the strips. Put in other words, this field 
describes what is left from the spin degrees of freedom of the virtual gluons forming the 
dense planar graph. 

It is convenient to parametrize the field which has 2(D — 2) independent compo- 
nents, as an antisymmetric tensor p^uiC) satisfying a number of constraints 

Vi±vVi±v = const, D = 2,3; 

( 6 - 5 ) 

ZVZV = const, p^ v p^ v = 0, D = 4 



PiivP^v = const, D = 2, 3; 

(6.6) 

PnvPnv = const, p^p^ = 0, D = 4 



Following the same logic as in the Polyakov quantization of the bosonic string |22 



we can 



look for a local action which is gaussian with respect to the x-field and such that on the 
classical trajectories the "intrinsic" normal tensor p„„ coincides with the normal element 
t\xv ~ Cabd a x^dbXv defined by the embedding. Since there is no intrinsic metric, the only 
possible quadratic action is given by the antisymmetric combinationi J d 2 ^p^ v {^)t^ v {^). 
The integration measure over the Grassmann manifold will be represented as a flat mea- 
sure for the antisymmetric tensor field with the constraints (|6.6| ) imposed through 
Lagrange multipliers p(£) and #(£)• In -D = 4 dimensions the Euclidean functional integral 
in the space of surfaces reads 

VS = Vx tl Vp liV VpV0 e~ A (6.7) 

A = J cP£[ie ab p> tv d a Xpd b x v + ^p(p» v p» v + M 2 ) + idp^p^] (6.8) 

where the constant M with dimension of l/[area] is the string tension parameter. This 
functional integral is formally quite similar to the path integral for the Eguchi string |2~4" . 
Here the Lagrange multiplier fields restore the reparametrization invariance which was not 
present in the Eguchi model. Depending on the dimension, it is possible to add topological 
terms but we are not going to discuss this point here. 



6 The same form of the action was recently suggested by A. Polyakov in ref. p3fl , where the 
antisymmetric field B tlu (x) was considered as an external field defined in the embedding space. 
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Returning to the case D = 2 we find, after integrating over the field j> Mi ,(£) = 
e /j.vP{£)i the following functional integral which strongly reminds the path integral for 
one-dimensional particle: 



Z = J VxVpVpexp ( - J d 2 i[it ab e^pd a x^d b x v + ^p(0(p(^) 2 + M< 
Vx^dpexp 



(6.9) 



_1 J d 2^ (dix 1 d 2 x 2 - d 1 x 2 d 2 x 1 ) 2 | ^ 2 

This analogy becomes more evident in the gauge x 2 = £2 i n which the integral (|6.9| ) is 
gaussian. The classical equations 

tabdapdbX^ = 0, p = it ab t^ v d a x^d b x v / p = ±iM (6.10) 

allow discontinuities of the p-field which correspond to embeddings having branch points 
and folds. Whether to take into account or not the configurations with folds depends on 
the regularization of the path integral (|6.9|). 

The loop correlators defined with such a measure are Stokes functionals in the sense 
of ref.0, i.e., they have a nonsingular area derivative, which property is shared by the 
(nonrenormalized) Wilson loops in the gauge theory. 



7. Discussion 

We formulated a model of random surfaces (lattice strings) satisfying the backtracking 
condition ( |1.3| ) and the lattice analogue of the loop equation ( |1.6[ ). This model describes 
the strong coupling phase of the U(N) gauge theory in the limit N — > 00 as well as the 
1/N corrections. 

In the simplest case of two dimensions this lattice string model is equivalent to the 
gauge theory up to the continuum limit. Nevertheless this is still not a continuum string 
theory with a local action on the world sheet. The most important step towards the 
continuum theory, the translation of the random surface Ansatz from the language of 
cellular complexes to the language of differential forms, is still to be done. The main 
problem to solveis to reformulate the contributions of the puncture operators in terms of 
local fields defined on the world sheet. In other words, we have to look for a continuum 
version of the puncture operators ( |3.31|) and ( |3.32| ) defined on the 1-cells and 2-cells of the 
lattice using the language of differential forms. 

The weights (|5.5|) of the branch points consist of a dimensionless "topological" term 
and a term proportional to the infinitesimal area A. (The higher powers in A can be 
neglected in the continuum limit.) The term proportional to A is easy to interpret in 
the continuum limit: it corresponds to a point-like singularity of the intrinsic curvature 
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weighed by a factor (—1). The factor A comes from the measure over surfaces, as we have 
noticed in |L0| . 



The real problem consists in the interpretation of the dimensionless weights. They 
seem to be associated with puncture operators which are pure derivatives. These are the 
only singularities in the "topological" limit A = of the model. In this limit the vector 
potential is restricted to be a pure gauge, A fl (x) = U~ 1 (x)d fl U(x). In this "topological" 
string theory the sum over surfaces bounded by a loop is always 1. 

Presumably there can exist only a finite number of local excitations with the lowest 
dimension. We have to study the classes of universality of local operators on the world 
sheet to extract what is left from the puncture operators in the continuum imit. 
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Appendix A. Expansion of the potentials in terms of moments 

A.l. The potential for the H-field 

The function F can be determined using the Ward identity 

[tr(0/0Ht0/0H) - N]e N2F(u) = (A.l) 

which is the matrix analogue of (|2.12| ). In the large- N limit this function was calculated 
explicitly by Brezin and Gross |25|1 . in terms of the eigenvalues of the Hermitean matrix 
Ht H. The result is given by two different analytic expressions depending on the strength 
of the external field H. These correspond to two different saddle points of the integral 
( p,19| ). In the "strong coupling" regime of small field H characterized by the condition 
tr(HtH) -1 / 2 > 2N the solution is 

o N i 3 

F ( R ) = n S ^+~ c - 2iV^ E + a/aTT^) ~c- z (A.2) 

i=l i,j 

where Xi,i = 1,2, N are the eigenvalues of the Hermitean matrix HJH and the param- 
eter c is determined by 



tr 1 



N VHtH + ci N^^Tc 



- T -7=^= = 2 (A.3) 

AT L^i . \ . T ~ v > 
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In the "weak coupling" regime defined by the condition tr(H"i~H) _1 / 2 < 2N the solution 
is given by the same formula (|A.2|) with c = 0. 

For sufficiently small field H the potential F(H) can be expanded as a series in the 
momenta 

a n = |(HtH)», (A.4) 



F l<*\ = E E *S^P^ (A.5) 

n=l fc!,...,A; n >l 



as 



The coefficients of the series were found in the large iV limit by O'Brien and Zuber [26 
a solution of a system of algebraic relations suggested originally by Kazakov in || . These 
relations are equivalent to eq. ( |A.1| ) written in terms of loop variables 



(^na n _i<9 n + ^2 [(n + k + l)a n a k d k+n+1 + ^a n+fc _i<9 n <9 fc ] - l) e N2p[a] = (A.6) 

n>l fc,n>l 



where we denoted «o = l,9 n = d/da n . Eq. ( |A.6|) is sufficient to determine the 1/N 
expansion of the link- vertices f[k 1 ,...,k n ]- I n the large N limit the identity (|A.6|) is equivalent 
to a system of recurrence relations [B] 



/[i] = 1 



fci-i 



/[fci,...,fc»i + y^fej/rAiiT+fc^fca A?,- fcj + /[*i-*,L\a]/[M] - ° 

j=2 fc=l q£P(L) 

where L = k2, k n and P(L) is the set of all subsets of L, including the emply set, and 
the bar means omitting the argument below it. The equations for the weights /rw of the 
disks § 

n-l 

f[i] = l; /[n] + X] /[*]/[»-*] = °' w = 1, 2, ... (A.8) 
fe=i 

are solved by the Catalan numbers 

/w = (-) ( - 1) ^ 1 ^ (A-9) 
L J n!(n — 1)1 

A. 2. The potential for the V field 

The potential S\(V) is defined by the analytic continuation of the r.h.s. of (|3.6| ) from 
the hyperplane V = U, ut U = I to the space of all complex matrices V. For this purpose 
one has first to express the characters ( |3.7|) in terms of the moments of the matrices V 
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and vt by replacing 
a series in the moments 



uf — > vf , u i n — > vf. After that the r.h.s. of 



(3 n = ^V, P- 



iV 



Vl n ; n = l,2,... 



n=l fci,...,fc„/0 



can be written as 

(A.10) 
(A.ll) 



The coefficients in the expansion ( A.ll ) of the heat kernel action can be found order by 
order by taking the logarithm of this series 



ri 

= 1 + N(trV + trVt)e" A/2 + (N 2 - l)(trVtrVt - l)e" A + ... 
5 A (V) = e- A / 2 (/3! + P-i) - e-\(3 1 (3- 1 + N 2 ) 

+ + /3_ 2 )(coshA _ iVsinhA )e -> + I (/3 2 + /£ J(-MdnhA 
+ 2(iVsinhA)2 )e -A + ___ 
which gives, in the limit N —>■ oo, 



(A.12) 



(A.13) 



Hi] 



-A/2 



S[2] = (coshA _ iVsinhA) e - A = (1 - A)e _A + <D(± 



N 



(A.14) 



si,-i = -e 

S[1>1] = (-iVsinh-1 + 2(sinh--)^)e- A = (-A + ^)e" A + O(^) 



^ ■ - ^ n2\ -A / a , 1 X 2\ -A 1 



-jy ■ v 2 JV' v " ' 2" y " ' " V N 2 ' 

In the continuum limit A — > the expansion (|A.11|) takes asymptotically the form 

Sx=o(P) = E -(A* + - /? n /?-n) = J] - T ^— T ^ (A.15) 



71 = 1 



which tends to the S- function 5(V, I) when V approaches the hypersurface vtv = I. 

The coefficients of the expansion ( |A.11|) satisfy a system of loop equations equivalent 
to the heat kernel equation on the U(N) group manifold ( |3.11| ) projected into the space 
of nonrestricted complex matrices 

2N d\ + tr ( va / av " Vta/avtj e ^ a ^(v) = o ( A16 ) 



7 This formula resembles eq. (22) in the paper by M. Douglas |fl4}| . We tried to establish an 
exact correspondence between his formalism and ours, but haven't succeed in this 
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Namely, 

nk 

V(nk)\n + k\fj n f3 k (J n+k + 

n=£0 k,n=£0 



(2d/d\+Y,\n\Pnd n + Yl [d(nk)\n + k\(3 n (3 k dn+k + ^(3 n +kd n d k ]j e^xtfl =0 (A.17) 



where we used the shorthand notation d n = d/d/3 n . 

This identity allows in principle to determine the coefficients S[fc lv . .,&„] m the expansion 
of S in the loop variables, but it is not triangular and does not produce simple recursion 
relations. Such recursion relations will be written below for the set of all possible loop 
correlators 

W kl ,..., kn = N n - 2 (trU k \..trU k -) connected (A.18) 
The coefficients S[ kl ,...,k n ] are related to the elementary Wilson loop correlators ([A.18Q 

by 

s [kl ,...,k n ]=W kli ... ikn , n^2 (A.19) 

and, for n = 2 

s [fci,fca] = W ki,k 2 ~ S kl ,-k 2 (A.20) 



In particular, s n = W n , as it has been noticed in |T0] 



A. 3. The one-plaquette Wilson loops 
Let us denote by 



W x [x] = -^log I VU e W [N 2 S x (U) + Nj2ynt?U r 

= J2 N2 ~ n h E w kl ,..., kn x kl ...x kn 



(A.21) 



the generating function for the connected Wilson loop correlators (|A.18p . 

Applying the identity ( |A.16| ) and integrating by parts in ( |A.21| ), we find the loop 



equation 

r d / 
2 ~g\ + E \ n \ x ™ d n + \ n Xn k Xk dn+k 

n^O n,k^0 fj^ 22) 

+ 1^6(kn)(\n\ + \k\)x n+k d n d k + )] e ^ 2w/ N = 
which is in a sense dual to ( |A.17|) . This equation is a particular case of the Migdal- 



Makeenko equations derived for a general Wilson loop in 
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In the limit N — > oo the loop equation (|A.22|) yields a system of recursive equations 
for the connected correlators VF/ Cl .. j c „ 



d n 1 n 

( 2 d\ + Yl k i) W kx,k 2 ,...,k n + J^J2 /, '' / U A •/, ,./,, A, k n 

i=l i=£j 

n ki — 1 

+ S \ k i\ W ki-k,L\c* W k,a = 

i=l fc=l aEP(L) 

In particular, for the simple Wilson loops W n = (l/N)(trU), eq. (|A.23| ) reads 

d 



2 — W n + nJ2w k W n _ k = 



k=i 



This equation has been studied in |fT9| . The solution of ( |A.24j) is 



< * l nm —I— I 1 1 V 



m=0 



(A.23) 



(A.24) 



(A.25) 
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Figure Captions 

Fig. 1. Three different ways of drowing the elementary surfaces associated with link- 
vertices. The edges are deformed for eye's convenience. 

Fig. 2. Elementary surfaces associated with plaquette- vertices 

Fig. 3. Two ways of drawing a piece of a surface obtained by identifying edges of surface 
elements. The link-vertex of order 2 creates a curvature 2n. 

Fig. 4. Vertices for the contact interactions between surfaces 

Fig. 5. A reducible surface 

Fig. 6. A piece of surface having a fold 

Fig. 7. Examples of Wilson loops: a) Circle, b) Flower, c) Pochhammer contour 

Fig. 8. A branched covering of the flower contour 

Fig. 9. Three ways of spanning the Pochhammer contour 

Fig. 10. A configuration of two loops 

Fig. 11. Microscopic structure of the world sheet in the weak coupling phase. 
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